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ABSTRACT 

The paper describes how, through simulation with physical and math- 
ematical models, the effects of ground winds on erected launch vehicles 
are being studied at the NASA, Langley Research Center. Experimental 
programs include both aeroelastically scaled models in wind tunnels and 
full-scale vehicles in natural winds. With the aid of a mathematical 
model having nonlinear aerodynamic damping, insight is gained which helps 
explain why in some instances the response at supercritical Reynolds 
numbers is a result of random forced excitation while in other instances 
is characterized by a self-excited periodic motion. 

Consideration is also given to the response of erected vehicles to 
atmospheric turbulence. It is shown that at the fundamental cantilever 
frequency of Saturn V class vehicles the gust velocities along the vehicle 
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are essentially uncorrelated for separation distances greater than a few 
vehicle diameters* Calculations indicate that the response is signific- 
antly reduced when these spatial correlations of gusts are 
account. 

; INTRODUCTION 

This paper will concern a search for new solutions to an old problem; 
namely, the effect of wind loads on flexible bluff bodies. For years civil 
engineers have grappled with the problem in the design of such structures 
as smoke stacks, suspension bridges, and tall masts* (See, for example, 
refs. 1 through U). I should point that included in these studies are 
contributions by the Director of this Conference, Professor Maher, who 
has conducted extensive wind-tunnel tests of suspension bridge section 
models. (See, for example, ref. 2 .) 

Starting with Goldman's investigation of ground-wind loads for the 
Vanguard in 1957 (ref.5), some form of wind-tunnel tests relating to 
this problem has been a part of the development program for practically 
every launch vehicle. Examples of such studies on specific vehicles 
are to be found in references 6 through 9, and other more general studies 
of vortex shedding from cylindrical bodies are given in references 10 
through 13. We find further evidence of the extent of effort presently 
being expended in this area by noting that at the annual AIM Conference 
on Structures and Materials held at Palm Springs, California in April of 
this year, five out of a total of seventeen papers In sessions on stru- 
ctural dynamics concerned the subject of vortex shedding from bodies of 
revolution. Despite this extensive research, there remain many important 
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questions to be answered, especially in the high Reynolds number range. 

With the aid of figure 1 let us identify the primary ingredients of 
the problem at hand. This schematic diagram shows a launch vehicle on 
a flexible support structure standing beside an umbilical tower. A 
steady wind imposes both static and dynamic loads on the vehicle. The 
static loads act primarily in the direction of the mean wind and the 
dynamic loads, which are associated with vortices shed from the vehicle, 
are largest in the direction perpendicular to the wind. These dynamic 
loads, which may result from either a random forced response or a per- 
iodic self-excited response, are frequently several times greater than 
the static drag loads. The adjacent tower structure and its turbulent 
wake may also have an influence on wind loads. These wind loads create 
problems in structural strength, guidance alinement, and clearance be- 
tween adjacent structures. 

The purpose of this paper will be to discuss recent studies in the 
area of ground-wind loads being conducted at the Langley Research Center. 
In particular, we will consider four items: wind-tunnel models and test 

techniques self-excited response, effects of atmospheric turbulence, and 


full-scale programs. 



WIND-TUNNEL MODELS MD TEST TECHNIQUES 


Aeroelastic Models 

The aerodynamic forces associated with ground-wind loads on erected 
launch vehicles are a result of flow separation from bluff bodies. Since 
such loads are critically dependent on Reynolds number, efforts are made 
to duplicate full-scale Reynolds numbers in the wind tunnel. Also, be- 
cause of the sensitivity of flow separation to minor surface protuber- 
ances, such as conduits or surface roughness conditions, careful attention 
must be given to the simulation of detailed geometric features. If steady- 
state loads were the only consideration, a rigid model having the same 
geometric features and Reynolds numbers as the prototype would be ade- 
quate for predicting ground-wind loads. Unfortunately, launch vehicles 
are not rigid and usually the predominant aerodynamic loads are not steady. 
Therefore, the structural dynamic properties of the vehicle, that is, 
mass, frequency, and damping, also become important considerations in 
model scaling. Experience has shown that the predominant dynamic response 
of launch vehicles to wind-induced loads is essentially that of the first 
mode alone. The significant additional nondimens ional parameters to be 
matched in the wind tunnel are then the reduced frequency based on the 
fundamental frequency of the vehicle rn^D/U, the structural damping, and a 
mass ratio involving generalized mass of the fundamental mode to the mass 
of air displaced by the vehicle. 

A model which simulates both aerodynamic and structural dynamic par- 
ameters is referred to as an "aeroelastic" model. Wind-tunnel tests of 
such models are believed to provide the most direct and realiable means 
of predicting the response of launch vehicles to a steady wind. (See, for 
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example, refs, ik and 15). 

Shown in figure 2 are some aeroelastic models of specific vehicles 
which have been used in ground-wind-load studies in the Langley 16-foot 
transonic dynamics tunnel during the past 4 years. Notice that in most 
instances the vehicle models are tested in the presence of simulated um- 
bilical towers or service structure. These towers are usually scaled 
only with regard to over-all geometry. The exception Is the erector tower 
for the Titan-Gemini which had scaled frequencies. The models - vehicle 
and tower structures - are installed on a remotely controlled turntable 
so as to permit response measurements to be taken from any wind direction. 

By using a freon test medium, which has a kinematic viscosity of 
about one-fifth that of air, Reynolds number simulation can be approximately 
achieved for all of the vehicles shown except for the Saturn V. For this 
model, which is 3 percent of full-scale size, the Reynolds number in the 
wind tunnel was one-third of the full-scale value. 

Even when Reynolds number is matched, however, there are differences 
in the flow in wind tunnels as compared with that in the natural atmosphere. 
In the wind tunnel the flow is approximately uniform and steady ; in the 
atmosphere near the ground the mean wind varies with height and is gusty. 

An approximate means of taking into account wind variations with height is. 
to relate the wind-tunnel Velocity to a calculated "equivalent” velocity 
which produces the same steady-state base-bending moment as does the actual 
wind profile. A reasonable approximation of natural wind profiles has been 
simulated in wind tunnels by use of a transverse grid of various sized rods . 
(See refs. l6 and 17- ) A pilot study is presently underway at Langley to 
determine the feasibility of applying these techniques In ground-wind load 


studies . 
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Modeling criteria for simulating atmospheric turbulence in wind tunnels 
have also been successfully applied in reference 18. These scaling laws, 
however, are not compatible with those used in the scaling of aeroelastic 
models. Some further consideration of the scaling of natural turbulence 
in wind tunnels will be given in a later section of this paper. 

Response Measurements 

The principal response measurements in these studies are the static 
and dynamic bending moments in two planes at the base of the model, and 
corresponding accelerations at the nose. Although time histories of the 
quantities are recorded, the most useful readout system employed in a 
two-axis oscilloscope and Camera. As illustrated in figure 3, the out- 
puts of X and Y bending-moment strain gages are fed to the X and 
Y axes of the scope and the sensitivities of the two channels are made 
equal. A time-exposure photograph of the oscillograph screen produces 
a roughly elliptical pattern which defines the envelope of the maximum 
bending-moment oscillations encountered during the data sampling period; 
typically, the model encounters between 2,000 and 3,000 cycles of oscil- 
lation ^during this period, which, in terms of the full-scale vehicle, is 
equivalent to wind exposure times of the order of 1 hour. The vector from 
the wind-off point to the center of the elliptical pattern represents the 
static moment , and the longest vector which can be drawn from the wind-off 
point to the tangent on the ellipse is the maximum resultant moment. Rote 
that information on the correlation between the two -moment time histories 
is displayed with this type of data presentation. It has been found that 
the probability of the lateral and dragwise bending moments reaching max- 
imum values at the same time is very low. 
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Damping Devices 

Structural damping has been found to he one of the key parameters 
which governs the susceptibility of a vehicle to wind-induced oscillations. 

It is very difficult, however to control damping in a model as accurately 
as other parameters. A viscous damper, shown in figure U , has been devel- 
oped which permits precise regulation of the damping in a model. The de- 
vice works on the principle of the classical Lanchester damper. It consists 
of a series of lead slugs (auxiliary masses) that are free to slide on con- 
cave trays inside of a cylinder filled with viscous oil. As the cylinder 
vibrates in a horizontal plane, the oil moves relative to the slugs, there- 
by causing energy to be dissipated. The damping can be varied by changing 
the number of slugs or the viscosity of the oil. Figure '5 shows the increase 
in damping in the 0.03-scale Saturn V model obtained with eight 72-gram slugs 
in the damper. 

Various damping devices have been used or proposed as means of allevi- 
ating wind- induced oscillations of smoke stacks. Usually the dampers are 
connected to the structure by means of guy wires (ref. 3 and k) . Since 
guy wire supports are often not feasible in launch vehicle applications, 
it is of interest to consider the use of auxiliary mass dampers, which re- 
quires no external connections , as a possible method for reducing wind- 
induced loads on full-scale vehicles. Such devices presumably could be 
attached to the vehicle during high wind conditions and then be removed 
prior to launch. 

Theoretical performance characteristics of an auxiliary mass damper 
are shown in figure 6 as a three-dimensional surface. It is assumed that 
an auxiliary mass m is connected to the vehicle structure through a 
spring and dashpot c, as indicated schematically in the figure. The 
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vehicle structure is represented mathematically as a single-degree-of free- 
dom oscillator having a generalized mass M and natural frequency 
The case shown is for m/M = 0.05* The vertical axis in the figure is the 
damping ratio £ associated with the vehicle mode; this was determined 
from the real part of one root of the characteristic equation for the 

coupled system. Plotted on the other axes are or /w, , the ratio of un- 

d 1 

coupled frequencies, and c/2mw^, a nondimens ional viscous damping par- 
ameter. This figure simply illustrates a well-known result from the 
theory of vibration absorbers (see, for example, ref. 19) that when the 
auxiliary mass is "tuned” to the natural frequency of the system to which 
it is attached, and c/2m.u^ is near an optimum value, large energy dissi- 
pation can be achieved. The previously discussed Lanchester type damper is 

represented by the curve for tu /w = 0. Note that this damper has optimum 

d 1 

performance at c/Smw^ 0.5. The wind-tunnel research applications where 
precise control of damping in aeroelastic models is desired, the fact that 

is relatively insensitive to variations in c/2mm 1 becomes an attractive 
feature; however, in full-scale applications, where weight of the auxiliary 
mass may be an important consideration, tuned dampers offer more than an 
order of magnitude improvement in performance. 

SELF-EXCITED RESPONSE 


Various and conflicting theories have been advanced regarding unsteady 
aerodynamic forces associated with vortex shedding from circular cylinders.. 
In the supercritical Reynolds number range - the range of primary interest 
for the present problem - experimental data indicate that in some instances 
these forces are random and insensitive to motion of the structure (Fung, 
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ref. 10) and in other instances periodic and strongly dependent on motion 
(Buel, et al. , ref. 12 and Den Hartog, ref. 3). This section of the paper 
will be devoted to consideration of these apparent differences. In par- 
ticular, a simple mathematical model, involving nonlinear aerodynamic damp- 
ing, is presented which has behavior consistent with both of the above- 
mentioned phenomena. 


Experimental Results 

In recent wind-tunnel tests at Langley involving Saturn I-B and Saturn 

V ground-wind load models violent self-excited oscillations were encountered. 

6 

These oscillations occurred at Reynolds numbers as high as 1 X 10 based on 
the model's maximum diameter. On both models the instability appeared at 
reduced frequencies, based on the maximum diameter, in the neighborhood of 
0.2 which happens to be the Strouhal number of the Karman vortex street at 
subcritical Reynolds number,. 1 Similar self-excited response phenomena have 
been observed on large, lightly damped, steel smokestacks and the solution 
has been to add damping (ref. 3) or attach aerodynamic spoilers to the 
structure (ref. 20). 

Now it might be argued that these high response conditions observed at 
a particular wind velocity could be in the nature of a forced resonant os- 
cillation rather than an instability. Data presented in figure J , however. 


It should be noted that for these models the peak response appears to be 
a result of vortex shedding from the lower stages at a Strouhal number of 0.2. 
Response characteristics were relatively insensitive to changes in nose shape 
or surface roughness conditions; for other vehicles configurations, nose shape 
and roughness may be the dominant factor affecting response. (See ref. 12.) 
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serve to disprove this argument. These data were obtained on the 0. 03- 
scale Saturn V model. The figure shows the variation of dynamic bending- 
moment response with velocity for three values of structural damping. First, 
note that for the smallest damping of C = 0.00 4 , the response peaks sharply 
at a reduced velocity of approximately .'It, 5-. The sample time history shown 
for the peak response point is a nearly constant amplitude sine wave at 
the fundamental model frequency. When the damping is increased from 0.004 
to 0.008 the sharp spike is eliminated and the response time history is 
characterized by a random-amplitude constant-frequency motion typical of the 
response of a lightly damped system to a random-forcing function. This com- 
plete change in character of the response time histories as damping increases 
suggests that at a certain critical velocity the structural damping of the 
system is a controlling factor which determines whether the dynamic response 
is a result of self-excited or externally forced motions. In either case it 
is apparent that structural damping of the model is an important parameter 
to he simulated. Furthermore, the commonly used scaling laws (ref. 21) based 
on the response of systems to random or periodic forcing functions, without 
consideration of motion-dependent aerodynamic forces, would be inapplicable 
for these eases of self-excited response. Similar results are shown in ref- 
erence 12. 


Conceptual Model of Self-Excited Response 
Let us now attempt to interpret the mechanism producing these observed 
dynamic instabilities. Scruton, in reference 22, has made extensive wind- 
tunnel studies of industrial structures such as smokestacks and towers, and 
has reported similar instabilities which he attributes to a nonlinear aero- 
dynamic damping that becomes negative at a particular wind velocity. Scruton* 
results were obtained at subcritical Reynolds numbers; however, his findings 
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appear to "be applicable to the present studies where Reynolds number varies 
from 1 to 6 million. 

A launch vehicle structure and the aerodynamic loads associated with 
vortex shedding are represented conceptually by the block diagram pictured 
in figure 8. Let Fg(t} be the aerodynamic forcing function which is in- 
dependent of vehicle motion and F^ft) to a nonlinear aerodynamic damping 
force which depends on the vehicle motion. In the equations of motion, also 
given in figure 8, the left-hand side describes dynamics of the structure in 
terms of a natural frequency and damping ratio and the right-hand side ex- 
presses the aerodynamic forces in terms of a lateral force coefficient 

C_(t) and damping coefficient C.. For the present purpose it will be 
Li x 

assumed that C T (t') is a random function of time. Depending on the re- 

If ' " ''' — 

lative magnitudes of the structural and the aerodynamic damping coefficients, 
the response x(t) of the structure in the direction perpendicular to the 
wind can take on either of the two forms Indicated in the figure. When the 
damping is predominantly structural, the response will be the characteristic 
random-amplitude const ant -frequency motion typical of a lightly damped system 
driven by a random forcing function. With aerodynamic damping present in a 
nonlinear form such as measured by Scruton, the possibility exists for self- 
excited motions indicated in the figure by the constant amplitude sinusoidal 
response. A representative plot of the variation of Ch with reduced fre- 
quency (see fig. 9) shows that over a narrow range of fD/U near 0.20, the 
aerodynamic damping has a destabilizing influence which depends strongly on 
the amplitude of motion. At this critical value of reduced frequency the 
destabilizing effect becomes smaller as the amplitude of motion increases. 
Thus , an energy balance is eventually established between the energy input 
by negative aerodynamic damping and the energy absorbed by positive 
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structural damping. This leads to a limit-cycle oscillation whose amplitude 
is inversely proportional to structural damping as shown in plot on right 
side of figure 9* 


Analog Computer Studies 

Further insight into the behavior of the nonlinear mathematical model 
presented in figures 8 and 9 may be gained from analog computer solutions. 
For this purpose the equation in figure 9 has been programmed on an analog 
computer in the manner shown in figure 10. The random lateral forcing 
function C L (t) was obtained by passing the output of a Gaussian white-noise 
generator through a first-order low-pass filter. The filter time constant 
was adjusted to make the power spectrum of C-j-^t) approximate closely the 
spectrum for fixed cylinders presented by Fung in figure 11 of reference 
10. Output of the noise generator was adjusted such that the root -mean- 
square of C T matched Fung's measured value of 0.12. As an approximation 
to Scruton' s results in reference 22, the nonlinear damping coefficient 
was assumed to be inversely proportional to a running average of the ab- 
solute value of response. 

The results are shown in figure 11 as a plot of maximum response 
observed during a computer run against damping of the structure. The plot 
on the left of the figure represents the system response with the motion- 
dependent aerodynamic force set equal to zero. This response varies in- 
versely as /if which is the relationship predicted for lightly damped 
systems acted on by a random forcing function. (See ref. 21.) According 
to Scruton's data this is representative of conditions when the natural 
frequency of the system is not in the vicinity of the Strouhal frequency. 
The right-hand side of the figure shows corresponding results when the 
natural frequency and the Strouhal frequency have approximately the same 
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values, say, f^D/U^O.2. Note the similarities in these analog time 
histories and the bending— moment time histories shown in figure 7 for the 
Saturn V model. 

The apparent inconsistencies between Fungus observation (ref. 10 ) - 
aerodynamic forces on cylinders are random and essentially independent of 
body motions - and the self-ex cited response observed for both wind-tunnel 
models and smokestacks might then be reconciled as follows: In Fungus in- 

vestigation, the cylinder motion was limited to reduce frequencies below 
0.12; Scruton (ref. 22) and the present wind-tunnel studies indicate that 
motion dependent forces, in the form of negative aerodynamic damping occur 
in a narrow range of reduced frequencies near the Strouhal number of 0.2. 

In order to shed further light on these and related questions , a generalized 
research study on two-dimensional cylinders will be conducted in the Langley 
transonic dynamics tunnel up to Reynolds number of 20 X 10. The Martin 
Company, George C. Marshall Space Flight Center, and Langley Research Center 
will be joint participants in this program. 

RESPONSE TO TURBULENCE 

Comparison of ’Wind-Tunnel and Atmospheric Turbulence 

Unlike the uniform profile of steady wind in a wind tunnel, atmospheric 
winds near the ground are typically as illustrated in figure 12. The pro- 
file patterns shown represent instantaneous distributions of horizontal wind 
that may occur as a field of turbulence is blown past an erected vehicle. The 
dashed line in the figure represents the mean wind averaged over, say, a 1- 
minute period; the little arrows represent more rapid wind fluctuations of 
wind components in the direction of the mean flow. As indicated in the figure, 
these unsteady wind components vary randomly in both time and space. Since 
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these winds produce random loadings on the structure, power spectral tech- 
niques appear to offer the most fruitful method of analysis and will he the 
approach discussed in the present paper. 

Before considering some mathematical techniques for predicting the re- 
sponse of launch vehicles to ground wind turbulence, let us first discuss 
the possibility of simulating such turbulence in a wind tunnel. Consider 
the question: To what extent would the turbulence structure representative 

of wind tunnels need to be modified in order to simulate "atmospheric" 
turbulence in wind-tunnel testing of aeroelastically scaled models? Figure 
13 will serve to answer this question. The figure shows typical plots of 
the power spectrum of wind-tunnel turbulence and atmospheric turbulence 
near the ground. The curves are plotted against the familiar frequency 
scaling parameter fD/U. Turbulence, whether it be in the earth's boundary 
layer or in a wind tunnel, can be characterized by a mean eddy size or 
"scale” length L, and the ms intensity of fluctuation components re- 
lative to the mean wind velocity, \f u^/U. For atmospheric winds near the 
ground, the scale of turbulence is of the order of the over-all length of 
Saturn V vehicle, or 10 vehicle diameters. In wind tunnels a typical value 
of L might be 1/10 the diameter of a ground-wind-loads model. Also, the 
intensity of turbulence is an order of magnitude different for the two cases, 
being, say, 20 percent of the mean wind for the atmosphere as compared with 
1 or 2 percent for wind tunnels. In summary, the figure indicates that over 
the frequency range of interest in the present problem, typical values of 
the scale and intensity of turbulence in wind tunnels and in the atmosphere 
are vastly different; therefore, study of the effects of gust loads assoc- 
iated with atmospheric turbulence in wind tunnels must await the development 
of new techniques for the generation and control of wind-tunnel turbulence. 


.4 
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Basic Power Spectral Relationships 
In absence of suitable wind-tunnel techniques for studying the re- 
sponse of vehicles to atmospheric turbulence, we must work either with 
mathematical models or actual structures exposed to natural winds. Both 
of these approaches are being pursued at NASA-Langley. This section of 
the paper will concern an analytical approach to the problem. 

Power spectral techniques have proved to be powerful tools for analy- 
zing the response of systems to random inputs. Such techniques for air- 
craft gust-load prediction have been under development for more than a de- 
cade and offer an excellent foundation on which to formulate the present 
problem. Of particular value is the paper by Houbolt , Steiner, and Pratt 
(ref. 23) which presents an exhaustive review and extension of power spec- 
tral techniques in relation to the response of airplanes to atmospheric 
turbulence. A basic equation for the power spectrum of a linear system 
acted upon by multiple random inputs Is presented in references 23 and 2k 
wherein the structure is assumed to he divided into an arbitrary number of 
segmented areas with a random forcing function acting at the center of each 
segment . 

* r U) = + « 22 h 2 h 2 * + * 33 h 3 h 3 * + ... 

+ 2H e (t 12 H i* H 2 + * 12 h\ + ••• 4 23 H 2* H 3 + ■ ■ • d (l > 

where 

$ (.«) power spectrum of a response variable r which may denote 

base bending moment, tip deflection, etc. 
cross spectrum of turbulence velocities at i th and points 


on the structure 



-xvi ii -i 6- 


IL frequency response function which gives the response due to 

a unit sinusoidal gust component velocity acting at the 
segment of surface area associated with the i** 1 point 

$ 

1L (to) = H^(-to) complex conjugate of H. 

Re denotes the real part 

In its general form this equation indicates that a considerable amount 
of information regarding the spectral description of turbulence is required 
in order to define the response spectrum. Specifically, the power spectrum 
of horizontal components of turbulence must be specified as a function of 
frequency at each of n points on the structure together with the cross 
spectra between all combinations of pairs of the velocity components. Since 
such information is rarely available in the form required, various simpli- 
fying assumptions must be made in order to obtain numberical solutions to 
equation (l). The most frequently made assumptions are that the gust ve- 
locities are in phase at every point on the structure and that the tur- 
bulence field is homogeneous; that Is, 


11 


12 


;= $ 


13 


In 


( 2 ) 


With these assumptions, equation (l) reduces to the following simple 

form 

♦ r (u) = !h 1 2 ® 11 (3) 

where 

H = + ... frequency response function for a sinusoidal gust un- 

formly distributed over the vehicle length 
power spectrum of gust components 

Bohne, in reference 25 , utilizes this equation in analyzing the response 
of launch vehicles to ground-wind turbulence. Various empirical expressions 
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for the power spectrum of turbulence, derived from data obtained from wind 
towers and airplanes, are available for use in such calculations. (See, for 
example , refs . 23, 26, 27, and 28..') 

An obvious shortcoming of the assumption of perfect correlation of 
gusts along the length is that dynamic response in the fundamental mode is 
likely to be grossly overestimated, especially if gust wave lengths corres- 
ponding to the fundamental structural frequency are small in comparison with 
the vehicle length. A more rational assumption would be to account for cor- 
relation of gust velocities on the basis that atmospheric turbulence is 
locally homogeneous and isotropic; that Is, its statistical properties in 
a given volume of air are the same at all points and are Independent of ro- 
tations of the reference axes. In an analysis of response of line-like 
structures to gusty winds, Davenport (ref. 29) makes similar assumptions 
and justifies them on the basis of experimental evidence. Under these 
conditions the cross spectra depend only on the separation distance between 
points and are independent of height above the ground. The following re- 
lations then apply 


$ = 
11 

*22 

*33 “ * 

.. $ 

nn 


$ = 
12 

*23 = 

*3^ = * 

n ,n+l 


$ := 
13 

*24 = 

$ = . 

35 

. . $ 

n ,n+2 

(4) 

and equation (l) becomes 



f $ 12 , * * 

L CH 1 H 2 + H 2 H 3 + • 


V“ )= *11 { H 1 H 1* + H 2 H 2* + 

. .. 2Re 



$ 13 * * 

♦Jf (H 1 H 3 + H 2 H !* 


- 1 * iJ H l\ + H 2 H n + l + 


)]} 


(5) 
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Cross Spectra of Turbulence 

In order to calculate the response spectrum on the basis of equation 
(5) the cross-spectra functions must first be specified* These functions 
are derived in appendix A from theoretical considerations which employ 
Taylor’s hypothesis (i.e. , fixed patterns of turbulence are transported 
at the mean wind velocity) together with the assumption that over the 
height intervals of interest the turbulence is homogeneous and isotropic. 
Results from the theory are plotted in figure 1 6 . Also shown in figure 
15(a) are experimental data obtained by Singer (ref. 30) from a U 00 -foot 
tower in wooded country at the Brookhaven Laboratory and by Davenport (ref. 

29) on a 500-foot mast in open grassland. It is encouraging to note that 
the experimental data from both of these sites are in reasonable agreement 
and that the theoretical curve falls within the scatter of the data. 

Probably the most significant feature of figure 15 is that for vertical 
separation distances greater than about 0.3 of U component wavelengths or 
0.5 of v component wavelengths the gusts are essentially unc or related. 

For example, in a 60-knot wind at a frequency corresponding to the Saturn 
V fundamental cantilever frequency (f 0.5' cps), correlation of u gust 
components is negligible for separation distances greater than two maximum- 
vehicle diameters. 

The ordinate in figure 15, defined as the square root of the coherency 
function, represents the modulus of the complex cross spectra which consists 
of a real (copower) part and an imaginary (quadpower) part. It should be 
pointed out that in the theoretical cross spectra the quadpower is zero where- 
as in the experimental cross spectra a small quadrature component was measured. 
These quadrature components are due in part to shear flow in the wind profile 
and also probably lags in the wind sensors. 





Applications 

In appendix B the frequency response functions indicated in equation 
(5) are derived in terms of natural vibration modes of a structure. The 
response spectrum for a two-dimens ional gust input Is then expressed as a 
product between the response spectrum for one-dimensional gusts and a so- 
called "attenuation" factor which accounts for the two-dimensionality of 
turbulence. For example, the spectrum of response in the x direction 
due to u gust components is (see eq. (B-15a) 


$ (k) 
x 


T (k)4> (k} 
un x 


, uni form 


where the subscript n refers to the mode number and T (k) is the atten- 

un 

uation factor given by equation (B-l6) in appendix B. Examination of this 
matrix equation indicates that T^tk] is dependent on such factors as ve- 
hicle geometry, mode shape, wind profile shape, and the cross spectra of 
turbulence between points along the vehicle. 

In order to present some specific results the gust attenuation factor 

T (k) = $ /$ . „ was evaluated for a Saturn V class vehicle in a 

un x x, uniform 

cantilevered unfueled condition. The Saturn V has a maximum diameter of 
33 feet and is about 350 feet tall. Natural frequencies of the first two 
cantilevered modes were assumed to be 0.48 cps and 1.85 cps. The steady 
wind profile shape considered followed a 1/5 power law variation with 
height as suggested in reference 31. 

Results of these calculations are presented in figure l£. Note that 
for reduced frequencies fD/U greater than about 0.1 attenuation due to two- 
dimensional effects is approximately twice as great for the first mode as it 
is for the second mode. A physical explanation of this result can be gained 
by recognizing that in the first mode a l80° phase shift for gusts on opposite 
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sides of the node line has a load amplifying effect. Similar results are 
noted in reference 2k. 

Although gust correlation effects tend to attenuate the first mode to 
a much greater extent than the second mode, the absolute response is, never- 
theless, predominantly that of the first mode. This result is associated 
with the rapid decrease in turbulent energy with increasing frequency. 

An interesting feature to be noted is that the curves in figure l6 are 
independent of the scale of turbulence . This is a consequence of the 
assumption made in appendix A that the gust wavelengths of interest are 
equal to or less than the scale of turbulence. 

Power spectra of first-mode response to u and v gust components are 
presented in figure IT for 1-D and 2-D turbulence. The response deflection 
has been made nondimensional by dividing it by the static deflection assoc- 
iated with the mean wind speed. In these calculations it has been assumed 
that the wind speed is the 99-9 percent probable maximum at Cape Kennedy 
(ref. 31) which is k3.k knots at the 2^0-foot height. The velocity at this 
height is the effective wind velocity which, if distributed uniformly 

over the vehicle, would produce the same static dragwise base bending moment 

as does the assumed wind profile. The mean square turbulence was assumed 
“2 —2 2 

to be u = v = 0.03 . And finally, the form of the power spectra for 

u and v was assumed to be given by equation (B— 1^4 1 with L - 10D Q = 

330 ft at all heights . 

Note in the figure that at the lower frequencies, corresponding to 
long wavelengths, the 1-D and 2-D response curves are almost identical. In 
the vicinity of the first-mode natural frequency, however, the effects of 
2-D turbulence are Indeed significant. From areas under the spectra curves 
the variance of response for 1-D gusts are approximately twice that of the 
2-D gusts. 
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FULL-SCALE STUDIES 


In previous sections of the paper we have discussed ways in which the 
response of launch vehicles to ground-wind loads can he simulated by means 
of models - both physical and mathematical. Whether or not these models 
are the intended analogs of their full-scale counterparts must he judged 
on the basis of comparisons between appropriate model and full-scale data. 
Although considerable effort has been devoted to wind-tunnel studies of 
aeroelastic ground-wind models , there is little in the way of full-scale 
data available for correlation with similar model data. 

Jupiter and Thor Vehicles 

In order to fill this gap the Langley Research Center is engaged in 
two programs aimed at obtaining ground-wind response data on full-scale 
vehicles. The first of these programs involves use of a Jupiter and a 
Thor vehicle which have been acquired as surplus property. These vehicles 
are to be used solely for ground-wind-load studies over a 1-year period 
at Wallops Island. Vehicle response , in the form of bending moments and 
accelerations , will be measured together with the wind inputs at two 
heights near each vehicle. In addition, other wind measurements will also 
be obtained from five elevations on a 250-foot tower located about l/H mile 
from the vehicle site. Special fast -response anemometers similar to the 
prototype shown in figure 18 (ref. 32} will be used. These instruments 
have flat response to frequencies several times greater than the fundamental 
vehicle frequencies. Vehicle response and wind input data will be auto- 
matically recorded on magnetic tape at regular intervals , or whenever the 
wind exceeds a preselected value. 
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The previously mentioned wind-tunnel tests of a 20-percent scale Jupiter 
model indicated that small st rake-type spoilers placed on the nose had a pro- 
nounced load alleviation effect. These results are shown in figure 19 as a 
plot of the maximum resultant bending moment against wind speed. One of 
our first objectives in the full-scale program, then, will be to determine 
if spoilers are in fact as effective as predicted in the wind-tunnel studies. 


Titan-Gemini 

The second full-scale ground wind response program has been established 
for the Titan-Gemini at Cape Kennedy. This program is being conducted by 
Mart in-Baltimore under contract with the Air Force and Manned Spacecraft 
Center. It will be closely followed by the Langley Research Center while, 
in addition to having conducted wind-tunnel studies of the vehicle and 
erector system, is providing the previously mentioned fast-response ane- 
momenters for the wind-measurement phase of the program. The bending- 
moment response of the vehicle will be read out from load cells which are 
a part of the launch pad hold-down structure. For each vehicle there will 
be two opportunities to obtain response data when the vehicle is not 
sheltered by the erector tower; these are during a fuel loading checkout 
and immediately prior to launch. 

Although the system is not yet fully operational, some bending-moment 
response data were obtained by the Martin Company just prior to launch of 
the first Titan-Gemini vehicle. These data are presented in figure 20 to- 
gether with wind-tunnel results and theoretical predictions of response 
due to turbulence. The only source of wind data available at the time of 
these measurements was a propeller-type anemometer mounted a few feet above 
the blockhouse and read out visually on a meter. Therefore, precise values 



for the mean wind and the intensity of turbulence at the vehicle are not 
known; the intensity of turbulence assumed in these calculations was 


Yu = V V^/U = 0.2. Although there appears to be reasonable agreement 
indicated between these full-scale measurements and the predictions based 
on wind-tunnel data and assumed turbulence conditions , quantitative com- 
parisons should not be attempted until the more precise measurements of 


the wind are available 
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CONCLUDING REMARKS 

This paper concerns the role simulation plays in prediction of effects 
of ground winds on space vehicles erected on the launch pad. It is shown 
that with present techniques, both wind-tunnel and mathematical models are 
required in order to predict the response of vehicles to gusty winds. By 
means of a conceptual representation , which involves a random forcing 
function and nonlinear aerodynamic damping, new insight is gained into the 
mechanism of self-excited response which has been observed at supercritical 
as well as subcritieal Reynolds numbers. In consideration of the response 
to ground-wind turbulence, it is indicated that the spatial correlation of 
gusts along the vehicle length can have significant effects. Finally, it 
is hoped that by correlation of response measurements now being obtained on 
full-scale vehicles with those predicted from models the adequacy of existing 
simulation techniques can be better established. 
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a ,b 
n n 


'L 

x 


c 

D 


D 


0 


f(.r), g(r) 


F s (t) 

F D (t) 


H (l) (i«) 
n 


K n (a) 


M M 
x, y 

P i, q i 


SYMBOLS 

generalized coordinates for the nth modal functions in 

the x and y directions, respectively 

lateral force coefficient 

aerodynamic damping coefficient 

viscous damping coefficient 

diameter 

vehicle maximum diameter 

nondimensional correlation functions of velocity com- 
ponents parallel and perpendicular, respectively, to a 
line r between two points in the field of turbulence 
frequency , ops 

random forcing function due to vortex shedding from a 
fixed cylinder 
aerodynamic damping force 

frequency response functions relating response to a unit 

sinusoidal input force at station i on the structure 

Bessel function of the third kind 

reduced frequency (k = idDq/Uq) , or spring constant 

Bessel function of the second kind for imaginary argument 

vehicle over-all length 

scale of turbulence (see eq. (A-8}) 

auxiliary damper mass or mass per unit length 

generalized mass of the nth natural mode 

bending moment about axis indicated by subscript 

orthogonal components of gust velocity at a point i 
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c,n %,n 


V * 1 


S nl ,S n2’ S n3 


t 9 r,T 


T ,T 
un s vn 


u s v 




U 

e 


W(z) 


x,y 

x , 
st 


z 

a = ^ 2 -n/^ 
Jb 

r 

<$ 

e 


generalized forces for the nth mode in x and y 

directions , respectively 

correlation distance 

real part of a complex quantity 

cross correlation of horizontal gust components at point 
i and j on a vertical line 
integrals defined in appendix B 
nondimens i on al time (s =■ U 0 t/D 0 ) 
time 

attenuation of response spectra in nth mode due to 

spatial correlation of u and v gusts (see eq. B-15 ) ) 

horizontal components of turbulence in direction parallel 

and perpendicular, respectively, to the mean wind 

amplitude of sinusoidal u and v gusts 

mean wind velocity, subscript 0 denotes reference Value 

equivalent uniform wind velocity 

function defined in equation (B-lt) 

deflection response or horizontal axes 

static deflection associated with mean wind 

deflections of nth natural mode of vehicle relative to the 

tip deflection 

axial distance along vehicle 

normalized cross spectra of gusts (see eqs. (A-15) and 

Ca-175) 

see equation (3-3) 

distance between segmented areas 
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C ah 
n n 


viscous damping ratio relative to critical 
damping 

s x n y n A 

Non dimensional response functions(^C n = n n - J 


rD„2- 


density ratio parameter (jc = ^ ® ^ 

n 

2 

constant in figure 9 (jc = 
nondimensional time delay (»■?) 

reduced frequency based on L 9 (^v = 
density of air or wind-tunnel test medium 


nondimensional separation distance 0 - 4 $ or 

standard deviation of a randan variable 




power spectra or cross spectra of a random variable 


circular frequency 


Notation: 

A dot (•) denotes differentiation with respect to time 
A prime (V) denotes differentiation with respect to s 
A bar (-) denotes average value 

An asterisk (*) denotes conjugate of a complex quantity 
Vertical bars | | denote modulus of a complex quantity 


^ j row matrix 


[] 

o 


square matrix 
column matrix 
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APPENDIX A 

DERIVATION OF CROSS SPECTRA OF VERTICALLY SEPARATED 
COMPONENTS OF TURBULENCE 


Consider the horizontal components of turbulence along a fixed vertical 
line representing the center line of an erected vehicle. In accordance with 
Taylor’s hypothesis it will be assumed that turbulence can be treated as a 
space pattern of velocities which are transported at a mean horizontal ve- 
locity U. In addition, the turbulence Is assumed to be homogeneous and 
locally isotropic. With these assumptions the correlation of velocities 
between two points depends only on the distance between the points so that 
a simple relationship exists between space correlations and time correlations 
The method of analysis to follow may be considered an extension of the 
cross spectra of vertical gust derived by Houbolt (ref. 2h) in connection 
with two-dimensional gust loads on aircraft. 

Cross Spectra of u(t) 

Consider the u component of turbulence at two points a distance Az 

apart on the vertical z axis (see sketch a}. 

Z 



Sketch a. 
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The quantity of interest is the cross correlation between u^(t + t) and 
u (t } from which the cross spectra can be derived. The cross correlation 
function, defined as 

R 12 (t] = lim | J u^CtJUgCt + x}at (A-l) 

T -*■ «> o 

can be expressed in terms of the one -dimensional correlation function of 
the turbulence field. The two such functions required are f(r) and g(r] 
whi ch denot e the c or re lat I on of velocity c omp on ent s parall el and t rans ve rs e , 
respectively, to a line r between two points in the field. (See sketch b. ) 



Consider the situation at time t. The velocity at point in sketch 
a is u^Ct ) . Upstream of point z^ a distance Ut is the velocity u^t 1 
which will arrive at t seconds later. These components, u^(t + t) 

and Ug(t), can be expressed in terms of components parallel and perpendic- 

^ r 2 

(Az) + (Ut) at time t. Thus 
u^(t + t) = p^(t } cos 6 + q^(t } sin 9 
u 2 (t) = p (t) cos 0 + <L 2 W Sin 9 


( A-2 ) 
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where 


sin 9 » 


A.z 


\J (Az } 2 + 


COS 0 .a 


UT 


s/ (Az } 2 + U 2 t 2 


Substitution of equation (A-2) into equation (A-l) yields 


RlgCt] = u 2 £f(r) cos 2 0 + g(r} sin 2 0^J 


(A-3) 


where f(r) is the correlation of velocities and in the direction of 

r and g(r) is the correlation of velocities q ^ and q^ in direction 
normal to r. These correlations are defined as 
Longitudinal : 


f(r) 


Rp 12 (r) p lPg 



(A-U) 


Transverse: 


g(r) 


(r) 


*12 


\l i0) 


h 1 


Also, in deriving equation (A-3) use has been made of the following relations 
which result from the assumption of isotropic turbulence 



and 

_ (A-5) 

p i q l = *2% = 0 

The cross spectrum of the u 1 u 2 components of turbulence is defined as 
the Fourier transforms of the cross -correlation function given by equation (A-3) 



% 
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12 


(w) 


r» t 

1 

* t. 


At 


(A-6) 


Thus in order to derive a specific analytical expression for the cross 
spectrum it is necessary to assume a specific form of f(r) and g(r) in 
equation (A-3). For this purpose the simple exponential correlation functions 
which closely follows the form, of isotropic turbulence observed in wind tunnels 
will be used. (See refs. 33 and 3^V) These functions are 


f (r) = e-W /T “ 


g(r ) = [ l--^]e 




:/L ^ 


J 


where L is the so-called Integral - scale of turbulence defined as 


(A-t) 


L = J f(r) ir 
0 


(A-6) 


Substitution of equations (A-?) into equation (A-3) gives 


R 10 (t) 


12 (.As) 2 + uV 

y (az) 2 +u 2 t 2 

[ U 2 x 2 e 


J (Az) 2 +U 2 t 2 


+ (Az) 2 L , /( + 


( A-9 ) 


Now let 


a - Az/L 
X « Ut/L 


Then 
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R 12 (X) = ^ C 1 " 




2J a + 




-Ja Z +\ Z 


{ A-10 ) 


and the cross spectra of turbulence then follows from equation (A-6) 


v W - ^ f (i - — ^ 

0 J - J7T?' 


2.2 . 

O +X — 10)T, 
e dx 


> 


( A-ll ) 


. (v)-^ r (1.-2! .) 

u 12 11 j V y 




2/d 2 + X 2 


e dX 




where v = wL/U. Integration of this equation produces the following ex- 
pression for cross spectrum 


-W o) ] 




"12 


-¥ ! C-!> 1 (u u.i-V 1 ’«‘«i] 


> (A-12) 


J 


where 


and 


K n (a) 

H n (l) (ia) 


a = a J 1 + v 2 


the Bessel function of the second kind for imaginary argument. 


the Bessel function of the third kind 
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Equation (A-12) is a real quantity indicating that the derived cross spectrum 
has a zero quadrature component . 

As the separation distance Az approaches zero, the cross spectrum (eq. 
(A-12)) reduces to the power spectrum of the u component of turbulence 

/T 

which, for f(r) = e~ ' , is simply 



(A- 13) 



Y u = l a [-2H 1 (1) CicO - aH^ 15 (ia)] (A-15) 


An important simplification results when it is recognized that the wave- 
lengths for frequencies of primary interest are in general less than the scale 
of turbulence; that is 2n(tJ},/w<L. The parameter a in equation (A~15) then 
reduces to a single nondimens i on al parameter 



rnAz 

U 


for 


U 


»1«Q 


Ca-163 


Therefore 
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S 2 - (i ^f)-^ H 0 (l > (i *§0 ] (A " 17) 

The absolute value of equation (A-1T ) is shown plotted in figure 15(a), to- 
gether with experimental data obtained by Singer in reference 30 and Davenport 
in reference 29. 


Cross Spectra of v(t ) 

The equation for the cross spectra of horizontal transverse v compon- 
ents of turbulence is identical to Houbolt T s expression for the vertical w 
components presented in reference 24 : . This expression is 



a ( 1 + 




a + 


V 2 ) 


3/2 L 1 


( 1 ) 


(ia)J | 


(A-1’8) 


In the limit for Az = a - 0, this equation reduces to the power spectrum for 
transverse components of turbulence which, for g(r) as given in equation 
(A-7), is 


2 2 
u 1 + 3v 


11 (l + v 2 ] 2 


CA-19] 


' V 12 


_ Y l^ i | -a 2 jjli/~4ig) ] - J } (A-20) 


11 


1 + 3v 


Therefore 
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Again, it can be stated that for the frequency range of interest in the 
present study (v >> 1.0) equation (A-20) is essentially a function of 
w Az/U alone and can be -written 


12 


TT m Az w (l) / . Azto ^ 

2 ~ H 1 V 1 —J 


(A-21) 


Equation (A-2l) is found to be an excellent approximation to equation (A-20) 
for wavelengths equal to or less than 2L. 

A plot of equation (A-2l) as a function of f Az/U is shown In figure 


15(b) , where f is frequency In cycles per second. 


In order to obtain somewhat simpler analytical expressions for use in 
response calculations , the cross spectra defined by equations .( A- 15 } and 
(A-2l) have been approximated by damped cosine functions. On the basis of 
a least- squares -fit , the following empirical equations were obtained* 


y = e 
U 12 


-H.A(f Az)/U ( f Az 


(¥*-) 


3 V U 


and 


-3.2(f Az)/U it ft Az > 

V 12 = 6 C0S 2\—J 


CA-221 
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APPENDIX B 


DERIVATION OF GOST RESPONSE EQUATIONS 


Frequency Response Functions 

Two sets of frequency response functions will be derived. One set re- 
lates the response of the vehicle in the x direction to a sinusoidal u(t) 
gust input and the other relates response in the y direction to a sinusoi- 
dal v(t ) gust input . The following assumptions are made: 

a. The vehicle is structurally and aerodynamic ally symmetrical about 
its z axis which is vertical. 

b. Wind forces on the vehicle are proportional to the local diameter 
such that the assumptions of two-dimensional strip theory are applicable. 

c. Wind forces are quasi-static (proportional to the instantaneous 
dynamic pressure) and act in the direction of flow relative to the vehicle. 

d. The unsteady components of wind are small in relation to the mean 
wind speed. 

e. Coupling between natural structural modes can be neglected. 

First, consider the deflection response of a vehicle in the direction of 
the mean wind. Let the x component of deflection be expressed in terms of 
natural modes of the structure 


x(z,t) = a 1 (t)x 1 (z) + a 2 (.t)x 2 (z) + ... 


(B-l) 


where x (z) is the nth natural mode shape and a (t) is a generalized 
n n 


coordinate determined from solution of the differential equation 

£ + 2, .A 2 a = V- 


CB-2) 



-XVIII -37- 


In equation (B-2) ^ is the equivalent viscous structural damping re- 
lative to critical damping, the natural frequency, the generalized 

mass, and Q x ^ the generalized force — each for nth mode of vibration. 

Assume that a gust velocity u(t) acts over an element of length e 
centered over station z^ on the structure. The generalized force in equa- 
tion (B-2) then becomes 

r l.,n = J Q l 2 P°D D [ U - Vn + ^dz CB-3) 

■where <$(z,z^) is a function which has the following properties 


«(z,z 1 ) = 1 when (z z - <_z ± (z 1 + J-) 


= 0 elsewhere 

With the assumption that the air density p is constant and the mean wind U 
is large compared with u and a Q , equation (B— 3) can be written 


i o r c- u-(z t) 

Q x,n = 2 pU 0 V L S ln - 2 t S 2n + 2 t ^ 


where 




nl 


[ w J-a(f) 

J n n U n 1 


S n2 = 


J w ( f> 
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With the notation 



9t 


J2. iL 

D o 3s 


equation (B-2) becomes 


p fi r u(z t] n 

5 n" + (2S n k n + 'V V + k n = K t 2 S m + T J CB-5) 


where k 


u> D. 
n 0 


3 ll> 


and a prime denotes , — . 

If the gust at z 1 is assumed to vary sinusoidally with time in equa- 
tion (B-5); that is , u n ^ z i s ^^ = u n o e ^ S 5 , then the frequency 

response function for dynamic deflections of a vehicle in the x direction 
associated with this gust input is 


H ,(k) 


!n P = 


K f-W n f Zl ] 


nO 


*h - k " + i( ^ 0 k n + K V k 


(8-6) 


Note that the static dragwise deflection has been subtracted from the total 
deflection in obtaining equation (B-6 ] . 

The response of the vehicle in the direction transverse to the wind can be 
similarly expressed in terms of natural modes 
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yfz 9 t) = b (t )y 1 U) + b 2 Ct)y 2 Cz} + 


(b-t) 


where b^ is determined from the equation 


+ u 2b = 


b + 2£ m b 
n n n n n n M 


(B-8) 


and y n (z) = x^ (z ) by symmetry. The generalized force for transverse re- 
sponse is 


r 

Ssn ” J 2 pG D D + l-”hi y n + v ( t ) 6 ( z » z 1 ) ] } y a sin Bdz 


Cb-9) 


where 


sin 8 


-by n + v6(.z,z 1 ) 


/ U 2 + [-b^ + v6(z,z 1 ) 


In a manner similar to that followed in obtaining equation (B-6) the 
frequency response function for the transverse response becomes 


v u 


nO 


•5 k f- W (z ) 
n l n 1 


v nO o 2 Z' 

UT k n +i C 2 S k n + 0 - 5 KS n 2 > ^ 


where n = . 

n D o 


For cases wherein the gust input is assumed to be perfectly correlated along 
the length of a vehicle the frequency response functions may be obtained by 
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integrating equations C.B-6). and ( B— 10 ) over the vehicle length. Thus for a 
uniformly distributed u(t) gust the frequency response function for de- 
flection in the x direction is from equation (B-6) 


H (k) 
^n 


x S 0 
n n3 


uniform n 


k + i(2?A + KS n2 )k 


(B-ll) 


and similarly for a uniformly distributed v(t) gust in equation CB-IO) 


H n fc) 


0,5k S 0 
n n3 


uniform k - k + i(2c k + 0.5 kS )k 
n n n ud 


(B-12) 


where S 


= f W (z)d(z/l). 

Ho J a 
3 0 


Response Spectrum 

With the turbulence spectra relationships in appendix A and the frequency 
response functions derived in the previous section, working equations, for the 
structural response spectrum for the nth mode can he formulated as follows; 

For a uniform one -dimensional gust input equations (.3) » (B-ll), and 
(B-12) combine to yield 
one-dimensional u gust 


2 

k S o2 
n n3 


uniform 


<V" - 


**“ 


(B-13a). 


one -dimensional v gust 
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uniform 


2 2 
0*25^ s* 

n n3 


*„/TT 00 


2 2 2 * 2 V/U ° 
K - k) + (2? A + °* 5K B S n2 ) k 


(B-13b) 


The turbulence spectra in these equations are expressed as a function of 
ojDq 

k = y- — and u and v are normalized with respect to Uq. With, these 
changes in variables the spectra given by equation (A-135 and equation (A~19) 
become > respectively 


u/XL 


(k) - 


2 A L A 

* U Q 2 1 + (L/D 0 lSt S 


(B-llta) 


4 v/U 


0 


(k) 


2 V 


L 1 + 3{L/D 0 ) 2 k 2 
D °[l + (L/D Q ) 2 k 2 ] 2 


(B-lUb) 


For two-dimensional gust inputs the following expressions can be derived 
two-dimensional u gust 


<£> (k) * T (k}* r (k) 

£ un £ n, uniform 

n * 


(B-15a) 


two-dimensional v gust 


<k) 


T (k)4 n (k) 
vn n n. 


uniform 


CB-15b) 


where the functions T (k) and T (k] are in the nature of "attenuation" 
un vn 

factors which account for the two-dimensionality of the gust field. These 
" attenuation" factors can be expressed in matrix notation as follows: 
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T n (k) = (£/ ^ 2 |_W n (j E ) j [y.^W ] {w n (je) } (B- 16 ) 

S nl 

where 

i = 1 , 2, ... m 
j = 1, 2, ... m 

The elements y. ,(k) , representing the real part of the normalized cross 
3- J 

spectrum between points i and j , are given for the u and v compon- 
ents of turbulence by the approximate formulas in equation (A- 22 1 . This 
function depends on the properties of turbulence and is independent of 
structural parameters. The quantity W( j e } , on the other hand, depends on 
the diameter and mode shape of the structure as well as the shape of the mean 
wind profile. (See eq. (B-UJ.) 
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Figure 1.- Factors contributing to ground-wind loads 
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Figure 3«- Time exposure of bending moments on an oscilloscope screen. 
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Figure 7.- Effect of damping on maximum dynamic bending moment response of 0-05 -scale Saturn V 

model. 
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Figure 8. - A conceptual representation of dynamic response of structures to wind-induced loads . 
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Figure 10.- Analog computer simulation of vortex- shedding phenomenon. 
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Figure 14.- Response equations for systems excited by multiple -random inputs. 
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Figure l6y- Effect of gust correlation on modal response of a launch vehicle. 
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.* Figure 17.- Power spectrum of vehicle response to atmospheric turbulence 
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Figure 18.- Fast response anemometer. 



Figure 19.- Effect of nose spoilers on bending moment response of 0.20 -scale Jupiter model. 
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Figure 20.- Response of Titan -Gemini to ground-wind loads. 



